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Abstract 

The dynamics of a particle which is linearly coupled to the environment of a Boson 
field is investigated. For velocity coupling and for local (position) coupling the reduced 
dynamics of the particle induces superselection rules for the momentum, if the Boson 
field is infrared divergent. Thereby the Hamiltonian of the total system remains 
bounded from below. 

1 Introduction 

Decoherence and superselection rules are the basis to understand how classical physics can 
emerge within quantum theory. The mathematical structure of quantum mechanics and of 
quantum field theory provides only a few superselection rules, the most important being 
the charge superselection rule related to the gauge invariance of the electromagnetic field, 
see e.g. [^] and the references therein. But there are definitively not enough of these 
superselection rules to understand the classical appearance of the world within quantum 
theory. A possible solution of this problem is the emergence of superselection rules due to 
decoherence caused by the interaction with the environment |^]. 

We investigate the dynamics of a particle which is linearly coupled to a Boson field. 
The Hamiltonian of the total system is bounded from below. This talk concentrates on 
a specific aspect of decoherence: decoherence due to infrared divergence of the Boson 
field. In the case of velocity coupling the infrared divergence induces a superselection 
rule for the momentum with uniform estimates for the emergence of the superselection 



sectors 1 12 ] . For the model with local coupling the infrared divergence also causes a strong 
decoherence between states of different momenta, but the superselection sectors are not 
uniformly induced. In this talk only some indications are given about how these results 
can be derived. The detailed calculations will be published in [13|. 

The models considered here have quadratic Hamiltonians, and they can be solved 
in an explicit way. Despite of their unrealistic simplicity such Hamiltonians have been 
successfully used to construct heat bath models and for understanding the approach 
to thermal equilibrium Q [Q], or for calculating decoherence effects due to radiation 
The structure of quadratic Hamiltonians is therefore rich enough for obtaining interesting 
results. 



2 The models 

We consider a massive spinless particle coupled to a Boson field of mass zero. The method 
presented here can be applied to a particle in with d > I. But to simplify the notations 
we restrict ourselves to one dimension. The Hilbert space of the theory is the tensor space 
(g) = >C2(M) (g) J^iTio) of the Hilbert space for the particle Tis = ^^(M) and of the 

^Talk presented at the Second Winter Institute on Foundations of Quantum Theory and Quantum 
Optics. S N Bose National Centre of Basic Sciences, Kolkata, India, January 2 - 11, 2002. 
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Fock space Hf = ^{'Ho) for the Boson field, where Ho is the one-particle Hilbert space of 
the Bosons. The Hamiltonian without interaction between the particle and the field, 

H(^o) = 1{P^ + ^^Q^)'^If + Is(^Hf (1) 

is the sum of the Hamiltonian of a harmonic oscillator or free particle (a; = 0) on the 
Hilbert space and the Hamiltonian Hf of a free field on the Fock space J^{Ho). As 

interaction we take 

- a velocity coupling Hp = P without harmonic force (a; = 0), 

- a local interaction Hq = Q (g) together with a harmonic force (a; > 0), 

where is the Boson field operator smeared with a suitably chosen test function h. 

For these quadratic Hamiltonians the field equations are linear and can be solved explicitly. 
To investigate the reduced dynamics we use the Heisenberg picture and take the algebra 
generated by the Weyl operators as algebra of observables. 

3 Weyl operators and classical phase space 

3.1 Quantum mechanics of the particle 

As Hilbert space of the particle we take Hs = >C^(K). The Weyl operators are defined as 
exponentials of the canonical variables, the position operator Q and the momentum oper- 
ator P, as Ws{a, b) = expz {aP + hQ) with the parameters (a, 6) G M^. In the momentum 
representation on the Hilbert space the Weyl operators are given by 

{Ws{a, h)f) [k) = exp l^-ab + iak^ f{k - b). (2) 

From this identity follows that the arguments (a, b) G M? can be interpreted as the position 
and the momentum variables of the classical phase space of the particle. The closure of the 
linear span of Weyl operators with the operator norm is the Weyl algebra. This algebra is 
smaller than the algebra of all bounded operators on Zl^(M). But the closure of the Weyl 
algebra with respect to the weak or the strong operator topology includes all bounded 
operators. 

For a one-dimensional harmonic oscillator with Hamiltonian H = ^ (P^ + uj^Q^^ , 
u; > 0, and unitary group U {t) = ex.p(—iHt) the time dependence of the canonical variables 
is U+{t)PU{t) = Pcosut - Qojsmut and U+{t)QU{t) = Puj-^smut + Qcosut. The 
dynamics of the Weyl operators then follows as U'^{t)Ws{a,b)U{t) = Ws{a{t),b{t)) with 

( :s ) - ^<') it)- - ( -rL. "itr ) 

where R{t) is a group of rotations on M?. For lo — > +0 the transformations R(t) have a well 
defined limit, which corresponds to the dynamics P{t) = P{0) = P and Q{t) = Q(0) + Pt 
or a{t) = a + bt and b{t) = 6 of a free particle. 

3.2 Field operators 

The basic Hilbert space for the Boson field is TIq = with inner product {f \ g) = 

J f{f^)9{k)d"k and norm ||/|| = {f \ /). The operator M is the one-particle Hamilton 
operator (M/) (k) = uj{k) f{k) with uj{k) = V^, k"^ = Y.l=i k = {ki, ...,kn) G M". 
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The operator M and its inverse are unbounded operators. To define real vectors we 
introduce the involution f*{k) := f{—k). The operator M is a real operator (Mf)* = Mf* 
with respect to this involution. 

The one-particle space TCq generates the Fock space J^{7iQ). Let us denote with fog 
the symmetric tensor product. Then the normalizations of the Fock space J-(Ti.o) are 
chosen such that the inner product of exponential vectors exp/ = l^^c + / + ^/ ° / + 
... G T{7io), f £ 7io, is (exp / | expg) = exp (/ | g). The normalized exponential vector 

exp (^f — ^ 11/11^) ) / ^ "^0) is usually called a coherent state. We use the notation cr{f) 
for the projection operator onto that state. 

The creation operators: a~^{h), h £ TLq, are uniquely defined by a'^ih) exp f := ho exp /, 
and the annihilation operators a{h), h S TIq are the adjoint operators a{h) = (a+(/i))^. 
These operators are normalized to [a(/), a^{g)\ = {f \ g)- 

The Hamilton operator Hp of the free field is uniquely determined by Hp exp / = Mf o 
exp/. The Hilbert space of real vectors f = f* £ Ti.o, is denoted by Wr. It is convenient 
to introduce a test function space <Sr C TC^ on which arbitrary powers of M are defined, 
"^R = {/ £ "Hr I M"/ G Tiu, n = —1, 0, 1, ...}. For real test functions / G <f^K we define the 
field operator ^{f) and the conjugate momentum field n(/) 



Hf) :=^K(M-5/)+a(M^5/) 



n(/) := z [H, cl>(/)] = ^ (a+(M2 /) - a(M2 /) 



(4) 



see e. g. |10|]. These fields are self-adjoint and they satisfy the canonical commutations 
relations [<^{f),^{g)] = [U{f),U{g)] = and [^{f),U{g)] = i {f \ g) . From these defini- 
tions follows that the arguments of the canonical fields n(ti) and ^{v) can be extended 
to vectors u G TL-m and v € TLm, where 'H±m are the real Hilbert spaces with the inner 
products (/ I g)^i = (/ I M^^gY The spaces Tim and Ti-m are dual spaces with respect 
to the inner product (/ | g) of T^k. 

The Hamiltonians M and Hp generate the unitary groups UQ{t) = exp{—iMt) on TCq 
and Up{t) = exp(—iHpt) on T{7iQ). The time evolution of the field ^{v) and of the 
canonical momentum field n(M) is 

^{v,t) = U^{t)<^{v)Up{t) = ^{cos{Mt)v) +U{M-^sm{Mt)v) , . 

n(n,t) = Uj{t)n{u)UF{t) = -$(Msin(Mt)u) + n(cos(Mt)n). ^ ' 

These identities are defined for test functions {u,v) G Ti—m x Tim- 
3.3 Weyl operators in QFT 

We define the unitary Weyl operators Wf{u,v), {u,v) G TC-m x Tim, by 

Wp{u, v) := exp {iU{u) + i<^{v)) . (6) 

The M-linear space of the arguments {u, v) G Ti-m x TCm can be interpreted as the classical 
phase space of the theory. The expectation of the Weyl operators between exponential 
vectors is 

(exp / I Wp(u,v)ex.p f) 
= exp(^\\ff -^\\u\\\-^\\v\\l + iIm(^M^u + iM-^2v\ f^j (7) 
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which includes the vacuum expectation. The expectation of Wf{u, v) in a K(ubo) M(artin) 
S(chwinger) state of inverse temperature /3 > is 

,r.r . ^^ / 1 / , , \ 1 / , , \ \ 

= exp 1^-- (^^/ I coth— - - l^v \ coth —v ^ j. (8) 

The dynamics of the Weyl operators follows immediately from as 

U;^{t)WF(,u,v)UF{t) =WF{u{t),v{t)) with 

^(*) ^ _ r,u^ ( ^\ r.u^ _ ( cos(Mt) M-i sin(Af t) \ 
v{t) J ^^>\v J'^^^l \-Msm{Mt) cos(Mt) J ^ ' 

The mapping R{t) is a rotation of the phase space Ti-mxTi-m with respect to the quadratic 
form (^u I M^uJ + | M~2v^ = + ||u||^ of the vacuum functional. 

4 Dynamics of the system particle and field 

4.1 The flow on the classical phase space 

The Hilbert space of the total system particle and field is Hs (^Hf = ® T{Tlo). 

The classical phase space can be represented by (a, u,b,v) S M x Ti^m x M x Tim and the 
Weyl operators are 

W{a,u,b,v) = expi {aP + bQ) If + Is (niu) + <^{v)) 
= Ws{a,b)(^WFiu,v) with {a, u,b,v) GRxn-mX-^x-Tim- 

The dynamics of the fields is uniquely determined by the dynamics of the Weyl operators 

U+{t)W{a, u, b, v)U{t) = W{a{t),u{t),b{t),v{t)). (11) 

where the flow {a{t),u{t), b{t),v{t)) on the classical phase space Mx'H_i]8 xMx'Hir is given 
by a one-parameter group R{t) = expLt of symplectic transformations. The generator L 
of the group R{t) is defined on the restricted phase space M x (Sir x M x <Sk and can be read 
off from the linear field equations of the theory. 

4.2 Reduced dynamics 

We assume that the initial state of the total system is given by 

P = PS^ PF eVinsC^np) (12) 

where ps G T>[7is) is the initial state of the particle and pF G T>[7iF) is the reference 
state of the field. The dynamics of any observable A € B{7is) of the particle is then 

A A{t) = $t [A] := tri.;7+(t) {A If) U{t)pF. (13) 

If A is the Weyl operator Ws{a, b) of the particle, then Ws{a, b) ® Ip = W{a, 0, 6, 0) 
is a Weyl operator of the total system as defined in Sect. 4.1 with 

(a, li = 0, 6, V = 0) G M X Ti-m. x M x Tim.. For quadratic Hamiltonians the time evolution 
C/+(i) {Ws{a, b) ® Ip) U{t) yields the Weyl operator (IllD W{a{t),u{t),b{t),v{t)). Thereby 



^The KMS states are temperature states on the Weyl algebra. They are defined also for Hamiltonians 
with a continuous spectrum, see [S. 
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the trajectories {a{t),u(t),b{t),v{t)) start from the initial conditions 

(a(0), n(0), b{0),v{0)) = {a, 0, b, 0). Since W{a, u, b, v) = Ws{a, b) Wf{u, v) the reduced 

dynamics (^) is simply calculated as 

[Ws{a, b)] = Ws{a{t), b{t)) x(a, b; t) (14) 

with the function xio-ib'it) = ivpWF{u{t),v{t))pp. A decrease of the trace x(O)^;0 for 
large t indicates decoherence effects due to the environment of the Bosons. 

If pf is the projection operator cj{g) onto the coherent state exp {g — \ , 9 G 'Hq, 

the trace x(o, 6;t) = tTFWF{u{t),v{t))a{g) is bounded by, see (0), 

|x(a,6;t)| = exp {--^ \\u{t)\\\ - i \\v{t)\\i^ . (15) 

The choice of the coherent states simplifies the calculations considerably. But the state- 
ments about decoherence are qualitatively not affected, if one takes another statistical 



operator as state of the environment |11] 



If the Boson field is in a state of inverse temperature /? > the function x(ffli&;i) is 
given by (|) 

Xp{t) = exp (-1 (u{t) I coth^n(t))_^ - \ (v{t) \ coth^t;(t)) J ^^^^ 
<exp(-iMt)||^-iHt)||?). 

Remark The technique to calculate the dynamics of a subsystem with the help of 
Weyl operators has been used by Davies to derive the approach of thermal equilibrium in 
a heat bath ^ . 

4.3 Evaluation of the dynamics and induced superselection rules 

4.3.1 Independent systems 

We first consider a free or harmonically bound particle on the real line and a free field 
with Hamiltonian (^). The operators 

Ro{t) = { , ] with Mo = f l) (17) 

y -Mosin(Mot) cos(Mot) j " \0 MJ ^ ^ 

are rotations on the phase space M x 7i-m x M x Tim. This rotation is the tensor product 
of (D and d). 

4.3.2 Interacting systems 

Velocity coupling The Hamiltonian with the velocity coupling 



H = ^P^ (g) If + P <E) <^>ih) +Is^Hf 
1{P0If + Is^ ^{h) f + Is {Hf - |$^(/i))^ 



is bounded from below if ||M~^/i|| < 1. The momentum P is a conserved quantity and it 



is possible to calculate the reduced dynamics of this model [12|. The Bosonic part of (18) 



is the van Hove model [|l^]. If h has low energy contributions such that h is no longer in 
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the domain of the operator M~2 ^ h ^ 'D{M~^), the ground state of the van Hove model 
disappears in the continuum, see |14] Q Sect. 6.1. Moreover, in the van Hove model and 
in our model the mean Boson number diverges due to the soft Bosons. As a consequence 
of this infrared catastrophe the selection rule of the momentum of the particle becomes a 



superselection rule |12]. Here we want to understand this effect with the help of the flow 
on the classical phase space. Assuming the stronger constraints 



M"^/i|| < 1, and h G P(M~t) c 



(19) 



the dynamics is unitarily equivalent to the free dynamics, and the flow can be easily 
calculated, see |l^. Starting from the initial values {a{0),u{0),b{0),v{0)) = (a, 0, 6, 0) we 
obtain 

/ ait) \ ( a + 6(M-i/i| {s\^Mt)M-'^h) ^ o?ht \ 



u{t) 
bit) 
\ vit) J 



bil 



V 



COS Mt)M-'^h 
b 

6(sin Mt)M-^h 

Since M has an absolutely continuous spectrum, the vector functions n(t) G 7i- 
vit) G remain bounded for t oo. 

These expressions have been derived under the conditions (1£). But the identities 
(pO|) remain meaningful under the weaker condition L/Vf~^/i < 1. If the low energy 



(20) 



and 



contributions in h are strong enough such that h ^ P(M~2)j there is a qualitative change 
in the time evolution: the Hamiltonian Jl^ ) is no longer unitarily equivalent to the free 
Hamiltonian, and the norms 



6^ II (/ - cos Mt)M-^h\\_^ = 6^ 
^ = 62 ||(sinMt)Af-i/i||^ = 62 



(/-cos Mt)M'2h 



(sin Mt)M~2h 



(21) 



diverge for t ^ oo b ^ 0. As a consequence the decoherence function decreases to zero 
in both the cases: an environment given by coherent states (^) and an environment given 
by a thermal state (|l6|). The formulas ( pO[ ) then imply the uniform bound 



[Wsia,b)]\\ <exp(-6V(t)) 



(22) 



with a positive non-decreasing function ifit) which diverges for t — > oo. If 6 = the 
relations (||) and (|o|) imply [Wsia, 6 = 0)] = Wsia, 6 = 0). 

Let jC?iR) be the Hilbert space of the particle in the momentum representation. For 
any interval I C ffi of the real line we denote with P(I) the projection operator, which 
maps into £^(1)^ i.e. /(x) e C^iM.) is mapped onto (P(I)/) (x) = fix) if x G I, and 

(P(I)/) (x) = if X ^ I. From (|) follows 



(23) 



P(l2)T^5(a,6)P(Ii) = if6^l2 



for any two non overlapping intervals Ii and I2 of the real line. The relations (|1J) ( p2D 
and (II) yield ||P(l2)$t [^^5(0, 6)] P(Ii)|| < exp (-(5V(t)) , where S = dist(Ii,l2) > is 
the distance between the intervals. Moreover any linear combination ^4 = ^ CjWsiaj, bj) 
of Weyl operators satisfies an estimate ||P(l2)<I>t [A] P(Ii)|| < C^exp [—6'^ipit)) with some 
constant Ca- These estimates imply: for all operators A of the Weyl algebra we have 



hm ||P(l2)«'t[A]P(Ii)|| =0, 



(24) 
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and for all bounded observables A the strong convergence 



lim P(l2)$t [A] P(Ii)/ = for all / G C'^{R) 



(25) 



follows. Hence the reduced dynamics leads to a superselection rule for the momentum. 

Remark A uniform estimate of the type (24) has been derived for all bounded ob- 
servables A with other methods in [IT2I. 



Position coupling For the Hamiltonian with the position coupling 

H = ^{P^ + uj^Q^) <S)If + Q<S) Hh) +Is®Hf (26) 



the flow on the phase space is - as in the free case - a rotation ( p^ R{t). But the square 
of the energy operator M is now given by 

M^- ^ ( ) . m 



The Hamiltonian (26) has a lower bound if and only if the operator (27) is positive. This 
condition is fulfilled if /i G Ti^ satisfies the norm bound ||M~^/i|| < uj. The operator (p7|) 
corresponds to the Hamiltonian of the Friedrichs model, see ||8|, and Sect. 6.2 in [Q]. 
If the function h satisfies some smoothness and support restrictions, the operator M has 
an absolutely continuous spectrum. Then one expects that either the harmonic oscillator 
decays into the ground state (if the Boson field is in a normal state), or the harmonic 
oscillator thermalizes into the canonical ensemble (if the Boson field is in a KMS state 
with positive temperature). But in the case ||M~^/i|| = (^h \ M^'^h^ = uP' the spectrum 
of the operator M includes zero. If, moreover the low energy contributions dominate, 
the mean Boson number diverges and strong decoherence effects emerge, which cause a 
non-uniform superselection rule for the momentum of the particle. The details of the 



calculations will be given in |13]. 
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